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September 25, 2002, which is hereby incorporated by reference for all purposes. 

FIELD OF TH F INVENTION 

The present invention generally relates to an apparatus and methods for determining 
absolute values of various properties of a physiological medium. In particular, the present 
invention relates to non-invasive optical systems and methods for determining absolute 
values of concentrations of oxygenated and deoxygenated hemoglobins (and/or their ratK>s) 
in the physiological medium. The present invention also relates to apparatus and methods for 
obtaining such absolute values by solving a generalized photon diffusion equation as well as 
its variations such as a modified Beer-Lambert equation. 

p Ar*rr,ROIJNP OF thf INVENTION 

Near-infrared spectroscopy has been used for non-invasive measurement of various 
physiological properties in animal and human subjects. The basic principle underlying the 
near-infrared spectroscopy is that physiological tissues include various highly-scattenng 
chromophores to the near-infrared waves with relatively low absorption. Many substances m 
a medium may interact or interfere with the near-infrared light waves propagating 
therethrough. Human tissues, e.g., include numerous chromophores such as oxygenated 
hemoglobin, deoxygenated hemoglobin, water, lipid, and cytochrome, where the 
hemoglobins are the dominant chromophores in the spectrum range of 700 nm to 900 nm. 
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According*, the near-infrared apectroacope has heen applied .0 measure oxygen levels in ,he 



concentrations. 

Various techniques have heen developed for the near-infrared spectroscopy, e.g„ 

wave spectroscopy (CWS). In a homogeneous and semi-infinite model, both of TRS and 
PMS have heen used to obtain spectra of an absorption coefficient and reduced scahenng 
coefficient of fire physiologica! medium by solving a photon dirmsion equation, and 
caieuiate concentrations of oxygenated ard deoxygenated hemoglobins aa weU as trssue 
10 oxygen saturation. CWS has generaUy been deaigned ,o solve a modified Beer-lamberi 
eq fion and ,0 meaaure changes in the concentrations of oxygenated and deoxygenated 



hemoglobins. u 

Despite their capability of providing the hemoglobin concensus as well aa the 
oxy gen samration, one major drawback of TRS and PMS is that the equipment ,s bulky and 
, 5 ex Lve. CWS may be manufactured a. a .ower cos, bu, limited in its utilriy because .. 

and deoxygenated hemog,obins. Accordingly, there is a neod for nove! CWS systems and 
as the oxygen saturation in the physiological medium. 
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SUMM ^PY ™ ™ ff INVENTION 

The present invention generally relates to an apparatus and method for obtaining 
absolute valuesofconcentra^^ 
their ratio, 

and methods for determining absolute values of oxygenated and/or deoxygenated 
hemoglobins in a physiological medium. . uAUr « 
in general, wave propagation or photon migration in a medium » descnbed by a 

generalized diffusion equation: 
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wave source and is a variab.e for an intensity of electtotnagnetic waves detected by a 
wavedetecto, Parameter is genera..y associated with the wave source and/or medmm 

Uaereof, mode of optica, conpUng between «te wave sonrce and medmm, and/or couphng 
daerebetween. Parameter «p» is generafly associate, with the wave detector and/or medtum 
and accounts for, e.g„ characteristics of the wave detector, optica. couphng mode betw^n 
.hewavedetecorandnredinnr.andtheassociatedconpbngios, Parameters-^ and p 
may a,so depend upon, ,o sonre extent, other system characters and optica, propen.es of 
.hentedinnt.inCndingnroseofcnrontopboresinc.ndennrerein. Paranteter -V maybe 
either a proportion* constant (including, e.g., 1.0) or a system parameter whtch may 

m edinm Parameter -B» generai.y aeconnts for .engths of optica, pate of photons or 

^eacwaves^oughthemed.nn.andisp— J bytheopbc, 
^ertiesof.eme.tnm.However.anexactva.neofparameter-B.mayasod^n 

su ch parameter "B» is convention* known as apath .engm factor. I. , appreetated ft* 
th epL«e r »B.maya 1S o^ t heva 1 neof..Owhere te generaUzed — et.nauon 

^LucedtotheBeer-Lamberte^on. To the contrary, parameter -L- tsgeneraUy 
leuy-dependent and accounts for a hnear distance hereon the wave source and wave 
detector. Parameter * may be either a proportional constant (inc.ndhrg, e.g., ..0) or 
system parameter which may be associated win, the wave sonrce, wave detector, and,or 
m edinm. Parameter «e, accounts for an optica, interaction or interference of photons or 
e.ectroma.ettc waves with an Mb chromophore inc.uded in the medium. ..is apprectat. 

Ifflcien, of me medium or the chromophores tnc.uded therein. VaHab.e »C, reptesents 
concentration of Ore i-ft chromophore inc.nded the medium, and parameter »o- ts etttter a 
proportionahty constant (inc.uding, e.g., 0.0) or a parameter which maybe assocated wdh 
L wave source, wave detector, and/or medium. Despite the numerous parameters of the 
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B eneralized diffusion equation (!) and various modified versions thereof which w», be 

ra,iOS,h rn f easpee t of, h epresen— .amedrodisprovrdedrosofveasefofwave 
^onsappulanop.ieafs.sfenr.avinsaffeasfonewavesoureeandarfe^ 

^ the physiological meditun ine.uding a, ,as, one chromophore, and 
W ave detector. Tb. wave equation, e. g ., the general diffirsron er.ua,, n (1), 

svs,emvariabies(e.g.,«f„»and-C i »)andsys tt mparame,ers(e.g., a, p, 

... .». » 6 » and-o"). Tbemeurodgeneraffyinofudesfiresfepsofobtongmntapfe 

s0UrC e-dependen« parameters (e.g., -a") and detector-dependent parameter (e.g P > 
n^ependen.andgeon.etiy-depe^enfparame^Ce.,, -B- and "L," reap— 
Cof intern!, e er.ua.ions, and obtaining absofufe values of ,be e—ns offb 

•P and » V )andfte medium- or ehromophore-dependentpaiametersCe.g., -a, ). 

This embodiment of the present invention offers several benefits over the pnor an. 
Contrary to the prior art CWS technofogy capab.e of measuring only the changes m the 

IfferentiL relative values of me physiological properties obtainabie by tire pnor «. CWS 
3„ PurtitenrtorcaswUlbedescHbedindetailbelow.theforegoingmethodofthe 
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Lb. of wave sources and detectors arranged in any arbitrary configurations. Therefore 
the foregoing method allows conation of optica, systems customized to sped* ctonca. 
applications without compromising their performance characlerrstics. 

Embodimentsofdusaspectofthepresent invention may include one or more of the 



following features. 



ing features. 

TT,e generalized difmsion equation (!) may be applied to an optical system wdh 
least one wave source and at least one wave detector: 



10 = 

where me subscript -m- and »n» present an m-tir wave source and an n,h wave detector, 

"^method may include me steps of app.ying equation (2) to the optical system to 
,5 obtain a firs, and a second set of equations, eliminating a. least one of a. P., * 5, - « *» 
the firs, and second se, of equations by performing mathematical operations .hereon ,o o 
afcds e t ofequatio„s,providm g a,,as,o„ecorre,ationwhichco— — ns 

iuc.udingB.and/orU,,— ng the above — inlo Ihe dmM ae. options 
20 ,0 replace such terms thereby, and obtaining an expression for absolute values of the 
concentrations of the chromophores (and/or ratios IhereoO based on known or measured 

Values Of Imnj Io,m? 8 i* . , A , 

The foregoing method may also inc.ude the steps of app.ying the opnca, system to the 
physiological medium including cells of organs, tissues, and body fluids, and meaaunng tire 
25 abide values of me chromophore concentrations (and/or their ratios) based on known or 
„> va,ues of U, and a,. The measuring slep may include an additional step of 
nrontoringconcentrations of oxy- or deoxy-hemoglobin, and/or a ratio thereof such as, eg., 

(tissue) oxygen saturation. f ^ mnr 
The foregoing memod may also include me step of determimng presence of tumor 
30 cells in a finite area of the medium or determining a presence of an ischemic cond.tion as 
W e.l in the alternative, the foregoing memod may also include the steps of applymg the 
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4 i oc r>f t 1 and &. For example, the method may be 

— — — 

and second set of electromagne p-exanote in the irradiating step, 

— — n^r^M.e.eseeondse.of 
ftefirstsetofelectromagnettcwavesmayhaveaflrstwave g 

elecroma^etic wave, may have a second wavelength which ts dtfferent from the 

one first ratio of two wave er.na.ions both of which are selected from one of the fits, and 

bn, to different wave detectors, thereby eliminating a„ Y, and o from the firs, ratio 

I lIX *e wave e q na.ions may be applied .0 two diffemn. wave sources bn. to .he 

Alternatively, me wave H The eliminating step 

, <, * tw of the first and second sets of the equations. A sum of or a 
^nmximating both parameters «y" and "5" as a unity. 

fonnnla of me medinm-dependen, and geometry-dependen, parameters as a po.yn^ 

geomefiy-dependen, parameters may be approximated as a constant 

to another aspect of me invention, an optical sys.em ts prov.de* to de.enn.ne 
a^nte vis of me eoncentmtions of ehromophores in the physical medmm and*r 
30 — Theoptiealays.emmayinclndeabody.asomeemodnle.adetectormodnle, 
Id mcessmgmodo.e. Th e aonme modnle is snpported by me body and , arranged ,0 
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electromagnetic waves having different wave charactensBc, The detector m 

e.echomagnetic waves transmitted through the tnedtum. T e 
5 arranged to operative* conp.e wi«t .he detector module and configured to solve a set of 

be designed to operate at a TRS.PMS or CWS mode. „,,„ riorart 
1 emhodimen, of ,he present invention offers severe, benefits over the nor ar, 

conventional TRS and PMS technologies, TRS and/or PMS opttca, systems of ^ 
systems of the present invention operating a. a CWS mode can — the abso.ute ahres 

device wearable by patients. ^ mnrPfl fthe 
Em bodimen,s of mis aspect of the present invention may mch.de one or more of me 



wavdength. Altemafive.y, me firs, setof e.ectromagnefic waves maybe «medby fi«* 
25 H wtve and me second - of e.ectromagnetic waves may be carried by a second earner 

different wave.engms, phase angles, amplitudes, harmonics, and men combmation 

of me chromophore concentrations (or .heir ratios) based on various variab.es and/or 
30 alters . L <* intensity of Cectromagnetic waves irradiated by the somce mod„.e, 
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parameters accounting for in.erae.ien or interference of photons or Ceotronragnene waves 

with the medium. . 

The wave equations may inc.ude a. leas, one term whieh is subs.anua.ly dependent 
„„ me optical properties of .he medium (i.e., medium-dependen.) and/or configuration of me 
; souree and de.ec«or modu.es (i.e., geome.ry-dependen.). Examples of such term may 
incta de, bu.no, limited*, »B- and «L- of equation (1) or -V and M.- ofequauon 
(2) The algorithm of me processing modu.e may include a. leas, one correlation expressmg 
a firs, function of me .erm as a second function of .he chromophore eoncen«mons (and/or 
rauos drereof). The second function may be any analytic function, e.g., a po.ynom.al of .he 
0 concentrations and/or ratios .hereof. Alternatively, the algorifhm may also be arranged fo 
approximate the second function as a constant. 

The source module may include at least one wave source and fhe defector module at 
.east «wo wave de.ec.ors. Alternatively, the source module may include at leas, two wave 
sources while fhe detector module may include a, leas, one wave detector. It - P-fcM 
, 5 however, mat both of .he source and de.ec.or modules include, respectively, a, leas, two 
wave sources and at least two wave detectors. 

m one aspect of medical application of the presen. invention, me foregomg ophcal 
sys ,ems and methods merefor may be used to measure .he absoluto values of — , ons 
of oxygena.ed and deoxygenatod hemoglobin and/or their ratio. Such opfca, systems wtl. 
20 beneficial in non-invasively diagnosing ischemic conditions and/or locating ischemta » 
various organ, and tissues. For example, .he optical system may be used to prognose or 
diagnose stroke, cardiac ischemia or other physiological abnonnalities originating from or 
charged by abnormally low concentration of oxyhemoglobin. Accordingly, presence 
of cancerous tumors may be easily de«ec.ed. The optical systems of .he presen. mvenuon 
25 may further be applied to cells disposed in epidermis, corium, and organs such as a lung, 
liver, and kidney. 

m another aspect of medical application of the presen. invenfon, the foregomg 
optica! systems and methods .hereto, may be applied to measure absolute values of .he 
concentrations of oxy- as well as deoxy-hemoglobins to diagnose vascular occ us.on durmg 
30 or after various surgical procedures including organ transplantation. In general, prognos, 
organ nansplantation depends on adequate supply of oxygenated blood to .ransplanted 
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organs during and post surgical procedure. The optica! system of the present invention may 
be used to detect vascular occlusion in transplanted heart, lung, liver, and tadney in ,ts 
earliest stage. 

In ye. another aspect of medical application of the present invention, the foregomg 
optica, systems and methods therefor may be applied to assess various absolute properties of 
,he physiological medium. Examples of such conditions may include, but no. tanned .o, 
concen.ra.ions (or .heir ratios) of lipids, cytochromes, water, and/or other chromophores m 
the medium. 

in another aspect of the invention, a method for determining an absolute value of 
concentrations of chromophores in a physiological medium is provided. A probe imda.es the 
physiological medium with electromagnetic waves a. a firs, and second wave.engms, where 
the firs, and second wavelengths are different The probe de.ec,s ,he efecfromagnenc waves 
emitted through me physiological medium. An absolute vafue of concentrations of 
chromophores is calculated torn .he de.ec.ed elec.romagne.ic waves. In some embod,men.s, 
.here are multiple probes a. different locations on the physiological medium for cahbration, 
generating regional information, monitoring changes over time (e.g., calculating tissue 

viability after surgery), and the like. 

The foregoing apparatus and methods of me invention may be employed for vanous 
applications, e.g., non-invasively disposed on the medium or, alternatively, to be mvaswely 
disposed on an internal medium. As used herein, me -chromophores" may mean any 
substiurces in a medium which exhibit a. least minimum interactioti with photons and/or 
electromagnetic waves emitting or propagating ureremrough. Examples of such 
cteomophores may include, but no. limited to, hemoglobins (e.g., deoxygenated or deoxy- 
hemoglobin (Hb) and oxygenated or oxyhemoglobin (HbO)), cytochromes, lip.ds, water, 
25 enzymes, hormones, transmitters, proteins, cholesterols, apoproteins, carbohydrates, 
cytosomes, blood cells, cytosols, and other optically interacting materials present m the 

animal or human cells. 

The terms -electromagnetic waves" include acoustic or sound waves, near-mftared 
rays infrared rays, visible lights, ultraviolet rays, lasers, and/or rays of photons. Unless 
30 otherwise defined, all technical and scientific terms used herein have the same meamng as 
commonly understood and/or nsed by one of ordinary skill in tire art to which this invention 
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belong, Although methods and materials similar or equivalent to those described herem can 
be applied and/or used in the practice of or testing the present invention, suitable methods 
and materials are described below. All publications, patent applications, patents, and other 
references mentioned herein are incorporated by reference in their entirety. In case of 
conflict, the present application, including definitions, will control. In addition, the 
materials, methods, and examples are illustrative only and not intended to be limiting. 

Other features and advantages of the invention will be apparent from the following 
detailed description, and from the claims. 

ttPTTTF nFSCRIPTI HN OF THE DRAWINGS 

FIG. 1 illustrates an example of a computer system that can be utilized in association 

with embodiments of the invention. 

FIG. 2 illustrates a system block diagram of the computer system of FIG. 1; 
FIG 3 shows an embodiment of an optical system of the invention. 
FIG. 4 is a schematic diagram of a sample optical system including two wave sources 
and four wave detectors according to the present invention. 

FIG. 5A is a plot of simulated values of G (i.e., a ratio of F, to F 2 ) at different 
wavelengths as a function of oxygen saturation according to the present invention. 

FIG. 5B is another plot of simulated values of G at different wavelengths as a 
20 function of oxygen saturation according to the present invention. 

FIG. 5C is yet another plot of simulated values of G at different wavelengths as a 
function of oxygen saturation according to the present invention. 

FIG. 6 is another plot of calculated oxygen concentration versus true oxygen 
saturation in a medium with a different background scattering coefficient and total 
25 hemoglobin concentration according to the present invention. 

FIG. 7 is a schematic diagram of a sample optical system including two wave sources 
and two wave detectors arranged in a linear fashion according to the present invention. 

FIG 8 is a time-course plot of total hemoglobin (HbT) concentration, oxygenated 
hemoglobin (HbO) concentration, and deoxygenated hemoglobin (Hb) concentration 
30 according to the present invention. 



15 
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FIG. 9 is a time-course plot of oxygen saturation according to the present invention, 
function of oxygen saturation. 

rr ,.„ MtsrRirn — urn mmnm 

The foHowing description provides preferred embodiments of non-invasive opttca. 
syst ems and methods for deterntinrag absohte va,ues of properties (or condtUons) of 
P Mo 1 ogica,tnedia I npa rt icu,ar,,hedescrip.io„prov i desopUcal systems and methods of 

L present invention for determining the abso.ute vahres of concentrauons of t 

hOT o g .obin concentration to total hemoglobin concentration which ,s a sum of to 

^ provides nove. methods of so.ving a Beer-Lambert equatton, a photon dtfiuston 
eauations and/or modified versions thereof. 

equation , exampleofacompotersystem « can be used ma— 

3 screens cabine.7 keyboard*, and mouse U. Mouse 11 can have one or more buttons 

ZZs incorporating computer codes that implement the invention, data for use wtth fiie 

irmemu.on.ercompn^re^.es.ragemeuiamc.udmgf.oppyd*^^ 
1 ry.syatemmemory.andharddHvescanbeuUH.ed. Additionany.ada.s.gna, 
I bod ^Laca ni er W ave(e,.,ina„e W otU„c 1 ud i ng,heh,teme,)canbe«hecom P u,er 

readable storage medium. „ nninl]tpr 
FIG. 2 shows a system block diagram of computer system 1. As m FIG. . , comput* 
system. mcludesmspUya.Iteyhoard^andmonsen. Computer system 1 ****** 
s bsystems such as a centra, processor 51, system memory 53, fixed storage 55 (e.g ard 

, Jera 63, and network interface 65. Other computer systems suttab.e for use wtth the 
i0 ^Lc^inCndeaddifionamrfewetsubsystems. F or cxamp.e, another computer 



PHOTP003 



11 



system could include more than one processor 51 (i.e., a multi-processor system) or a cache 
memory. 

The system bus architecture of computer system 1 is represented by arrows 67. 
However, these arrows are illustrative of any connection scheme serving to link the 

5 subsystems. For cample, a local bus could be utiHzed ,o connect processor 51 tomemory 
53 and display adapter 59. Computer system 1 shown in FIG. 2 is but an example of a 
computer system suitable for use with the invention. Other computer architectures havtng 
different configurations of subsystems can also be utilized. 

In one aspect of the invention, a nove. method is provided to solve the modified Beer- 

10 Lambert equation and/or tire photon diffusion equation applied to an optica, system includmg 
a source module and a detector module. The source and detector modules may tnclude, 
respectively, at feast one wave source and one wave detector. However, it is generally 
prefered fha. the source and detector modutes inctude at least two wave sources and two 

wave detectors, respectively. 
15 As described hereinabove, equation (1) is the generalized governing equation for 

describing migration of photons or propagation of electromagnetic waves in a medrum: 



20 



25 



/-.J-f/.-e' ' 



(1) 



Bis appreciated fha. the system parameters V and »5» may have the value of 1 .0 and »o» 
may be 0.0. One simplified version of equation (1) may be obtained when .he parameters 
AyA and A5A are approximated as a unity: 



The 



-m-IUgM (3a) 
/=a^•/ 0 •e , ' ' 

conventional "photon diffusion equation" has the same form as equation (3a): 



I=SDI„e [ ' 



30 



(3b) 
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where -S- corresponds .0 »«» of equation (3a) and generally accounts for charactensucs of 
Ure wave source such as power and configuration thereof, mode of optical coupfing between 
.he wave source and medium, and/or optica! coupting loss therebetween, -D- corresponds to 
... of equation (3a) and generally accounts for characteristics of the wave detector, mode of 
5 opticalcoup.mgbebweanmewavedetecwandmeutum.andyormeaasociatedcoupbng.oss, 

md «A» corresponds to »o» of equation (3a) which may be either a proportionality constant 
or a parameter associated with the wave source, wave detector, and/or medtum. 

For illustration purposes, an exemplary optical system may include, e.g., two wave 
sources (SI and S2) each emitting elecfromasrefic waves of wavelength X, and two wave 
,0 de.ec.ors(Dl and D2) arranged to detect a, least a portion of such electromagnetic ^waves. 
Appiymg the photon difiusion equation (3b) to each pair of fire wave sources and detectors of 
the exemplary optical system yields the following set of equations: 



15 l£ ID , = JJS,D,e l ' 



-«iii«(Z«« J|< 5) £ «» M+il 



20 



Is'lD! = IiS,D,e 

I 1 ( 4c > 

(4d) 



where me superscript X, denotes that various variables and parameters are obtained a. the 

25 wavelength of U- 

A simple mathematical operation may eliminate at least one system parameter from 
,he equations (4a) .0 (4d). For examp.e, the source coupling factors such as S, tntd S, may be 
canceled .herefrom by taking .he firs, ratio of .he equation (4a) to (4b) and by fakmg fire 
fourth ratio of the equation (4d) to (4c). Logarithms of me firs, and second ratios are tite„ 

30 .aKenfoyiefdwhafareconventionalfytermedas -optica. densities- (i.e., <wf •.«■"*» 
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aiogarithmof 4,/^ and OD* definai as a bgarita, of 4, 

these optical densities are generally insensitive ,o exact modes of optical coophng between 

the wave source and the physiological medium: 



10 
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1 ^S1D2 ^2 



(5a) 



2 ^52D1 1 



Other system parameters may also be eliminated through reformulating the above 
equations (5a) and (5b). For example, the terms including the detector coupling factors, D, 
and D 2 , may be canceled by adding equation (5a) to (5b): 

OD x > = OD\ + OD l i = F A '£ s?C, ( 6a) 

i 

where F A| = ^ 1D! L S1D! - BSntW fBfe D ,L s!D1 - B' io2 L sim ) (6b) 



As matufes, in the equation (6b), F K is primarily determined by confignrations of Are wave 
sourees and detectors (i.e., «L=s" which are predominant* geometry-dependent and whtch 
acconn, for distances between each pair of a wave sonrce and a wave detector) as wel. as the 
path length factors (i.e., »B-s» which are predominantly medium-dependent and wmch are 
determined by the optical properties of the physiological medium and/or electromagnet* 

waves). . 

Equations (6a) and (6b) may be applied to the physiological medium m order to 
25 obtain quantitative physio.ogica. information such as concentrations of the chromophores 
and/or .heir ratios. Numerous substances contained or suspended in the medium may be 
capable of interacting or interfering with photons or electromagnetic waves imptngtng or 
ptopagating therethrough. However, in many physiological media, hemoglobins such as 
deoxygenated and deoxy-hemoglobin (Hb) and oxygenated or oxyhemoglobin (HbO) are the 
chromophores of the most physical interests. Applying equations (6a) and (6b) to such 
physiological medium yields: 



30 
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(7a) 



where [Hb] and [HbO] respectively represent concentrations of Hb and HbO. 

By arranging the wave sources, SI and S2, or additional wave sources, e.g., S3 and 
S4 to irradiate a second set of electromagnetic waves having a wavelength X 2 which is 
different from the wavelength X u a companion equation of the equation (7a) is obtained as 
follows: 

(7b) 

OP* le^Q e%,[Hb]+ elo[HbO] 
F h 

Accordingly, mathematical expressions of two system variables [Hb] and [HbO] can be 
readily derived from two algebraic system equations (7a) and (7b) as follows: 

£ HbO fit ~ £ HbO ph (g a ) 



b Hb 



OD h 



OD h 



[Hb ° ] ~- *«i,-«M8P 



(8b) 



25 



where F h = (B* D2 L SW2 - B$ im L slDl )+ (B^ slDl - B^ m L slDI ) 

F h = {Bl] D1 L sw2 - B? im L Sim )HB* m L™ ~ B* D1 L S1D1 ) (8d) 



and 



Expressions of other physiological properties, indices or variables may also be 
derived from Ute above equations. For example, oxygen saturation (SO,) is a frequently used 
30 index for diagnosis of ischemic conditions and is generally defined as a ratio of concentrate 
of oxy-hemoglobin to total concentration of hemoglobins (i.e., [HbT] = [Hb W HbO]) 
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so 2 = 



[HbO] _ [HbO] 



[HbT] [Hb]+[HbO] 



(9a) 



incorporating equations (8a) and (8b) to equation (9a) yields a following formula for the 
oxygen saturation as a function of the extinction coefficients (s=s), optical densities (OD=s), 
and medium/geometry-dependent factors, and F^ : 

10 ,*°*LtL_,- 

1W on 1 ' F h m C9b) 

It is noted that the extinction coefficients of Hb and HbO at wavelengths U and X 2 
can be obtained experimentally or from the literature, and that the optical densities may be 
15 readily measured experimentally. However, estimation of F* and is not straight- 
forward because the path length factors including the terms F* and F^ may depend on 
specific types of the physiological medium as well as optical or energy characteristics of 
photons or electromagnetic waves. 

It is appreciated that the absolute value of [Hb], [HbO], and/or oxygen saturation may 
20 be obtained by estimating F^ , F^ or their ratio as a function of, e.g., [Hb], [HbO], and/or 
oxygen saturation. For example, it may be assumed that F^ , F^ or their ratio may only 
marginally depend on background optical properties and configurations of the wave sources 
and detectors. It is believed that these assumptions are fairly accurate in linear optical 
processes such as migration of photons or propagation of electromagnetic waves in the 

25 physiological media. 

Once the correlations of the ratio of F* toF^ with oxygen saturation is obtained for 
different physiological media by simply measuring the optical properties thereof, such 
correlations may be incorporated into equations (8a), (8b), and (9b), and the absolute values 
of [Hb], [HbO], and/or oxygen saturation may be obtained. In particular, a ratio of F* to 

30 F^ may be approximated, e.g., as a polynomial of oxygen saturation as follows: 
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.M^Lt a SOl-e* (10) 
e » OD h r„ ' 1 » v 



where coefficients of each term (i.e., (i.e., «, <«, «3 • • •) may be obtained by, e.g., 
theoretical derivation, semi-theoretical estimation or numerical method best-fitting 
experimental data obtained between the values of G and oxygen saturation. By incorporating 
the formula for G of equation (10) into equation (9b), the absolute value of the oxygen 
saturation maybe obtained from known values of the extinction coefficients (i.e., •«-•) 
and experimentally measured optical densities (i.e., «OD=s») as follows: 



G= —j~= Z a t SO{ - a 0 + a t S0 2 + a 2 SO\ * a 3 SO> + ... (1 D 

F 2 j-o 

Equation (11) is generally solved numerically, but an ana.ytical expression for the oxygen 
saturation may also be obtained when only a few firs, terms of fire polynomials are adopted 
s „as,oappreximateG„i.e..merafioof FMoF', . Other methods may also be apphed to 
approximate G. For examp.e, G may be estimated as a function of [Hb] and/or [HbO], 
although it is noted that the accuracy of this estimation may depend upon the one-to-one 
correspondence between G and [Hb] and/or [HbO]. Alternative.* G may be approximated 
as a constant as well. This approximation may be a reasonable assumption when F \ and 
F \ are relatively constant or tend to vary in proportion to each other according to different 
values of [Hb], [HbO], and/or oxygen saturation. In the alternative, the value of » W may 
be adjusted by manipulating geometric configuration of the wave sources and detectors so as 
to render G stay constant or vary in a pre-deteimined manner. 

Similarly, each of f\ and F \ maybe approximated as a function of [Hb], [HbO], 
and/or oxygen saturation, tn the aftemative, F* 1 and F* may also be assigned specific 
values which may best approximate the optica, system and/or the physiofogicaf medmm of 
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interest. By taking the simplest approach of approximating F* 1 and F h to be a unity, the 
absolute values of [Hb], [HbO], and oxygen saturation may be obtained as follows: 

[/#>] = t h _ E h £ h 

£ Hb £ HbO b Hb t> HbO 



s\OD* - e h J)D^ (12b) 



E m b HbO b Hb°HbO 



so 2 = 



OD h 



(12c) 



"Hb 



£ Hbo) /nn ; i • 



OD 



20 



It is noted that [Hb], [HbO], and/or oxygen saturation obtained from the equations 
(12a) to (12c) (and/or other approximation methods described hereinabove) may be less 
accurate than those obtained from equations (8a), (8b), and (9b). Nevertheless, as long as the 
foregoing assumptions hold valid, one-to-one correlations may be expected between the true 
values of [Hb], [HbO], and oxygen saturation and those obtained from approximating 
equations (12a) to (12c). Such correlations may be determined once the optical properties of 
the physiological medium are known. For example, extinction coefficients, absorption 
coefficients, and/or scattering coefficients of the physiological medium (or those of the 
chromophores) may be determined for [HbT] and oxygen saturation. With known optical 
properties, oxygen saturation maybe estimated at different levels of [HbT] through 
simulations of the diffusion equations and/or through experiments. Equations (12a) and 
(12b) may then be used to back-calculate [HbT], and a correction function can be calculated 
which correlates the calculated [HbT] with the true [HbT]. Similar or identical approach 
may be applied to calculate correction functions for [Hb] and/or [HbO] as well. It is noted 
that these methods may be applied to different physiological media (e.g., different human or 
animal subjects) to assess different optical properties and, therefore, to obtain different 
30 correction functions. 



25 
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It is appreciated that the foregoing methods are applicable to any optical system and 
physiological media where migration of photons or propagation of electromagnetic waves 
may be reasonably described by the generalized governing equation (1). It should be noted 
that the parameter eliminating step of the foregoing methods may be applicable regardless of 
5 thespecificnumericalvaluesassignedtotheparameters «y<< and «8'<. For example, y can be 
eliminated by taking ratios of equation (4a) to (4b) and equation (4d) to (4c), and 8 can be 
eliminated by taking the ratio of and F* 2 ■ m addition, the foregoing method may also 
be readily applicable to any modified versions of the governing equation (1) where the 
optical interaction or interference of the medium is described by absorption coefficient or 
10 scattering coefficients of the medium (or the chromophores included therein). For example, 
by assigning an adequate value and unit to parameter »y,» such modified equations can be 
converted into equations substantially similar or identical to the governing equation (1). 
Accordingly, it is manifest that the foregoing methods may be deemed universal methods for 
solving the generalized governing equation (1) for the chromophore concentrations and/or 
15 their ratios. 

It is further appreciated that the absolute values of the chromophore concentrations 
(or their ratios) may be obtained by variations of the foregoing methods. For example, the 
detector coupling factors, D, and D 2 , may first be eliminated from equations (4a) to (4d) by 
taking the third ratio of the equation (4a) to (4c) and the fourth ratio of the equation (4d) to 
20 (4b) as follows: 



4 h\m l S\ ^1 ZD 



(5d) 



30 



Similar to equations (5a) and (5b), this variational method yields optical densities. OD 3 and 
0&\ which are substantially insensitive the coupling mode between the wave detector and 
medium. By adding equation (5c) to (5d), the logarithmic ratios (i.e., one ratio of intensities 
of electromagnetic waves irradiated by the wave sources and the other ratio of the source 
ipling factors, S, and S 2 ) also cancel each other, yielding: 



coui _ 

19 
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5 



where K = i*& u L m - Bl\ m L stDl )HB? wl L swl - B? 2D1 L S2D1 ) (6d) 

By applying equations (6c) and (6d) to the physiological medium including oxy- and deoxy- 
hemoglobins, following equation (7c) is obtained: 

nn'< C^ c ) 



10 



F 1 

1 34 



Similarly, a companion equation of equation (7c) may be obtained by applying the second set 
of electromagnetic waves having a wave length X 2 : 

od^ ( 7d > 

r 3 4 



Accordingly, by solving equations (7c) and (7d), mathematical expressions of two system 
20 variables [Hb] and [HbO] can be obtained as follows: 



[Hb] = "i; ^ _ A2 p h 

t Hb t HbO °Hb°Hb 



25 _q^__ OD£ 



<4~^r (8f) 



30 



[HbO]= h h h i 

£ Hb b HbO C Hb C HbO 



where F » = " 5 -- L -> )+ (2? ™^"» " ( 8 §) 

and K 2 = ^\ m L sim - B? Wi L stm )+ (B>] m L SW2 - B? 1D1 L S2D2 ) (8h) 
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The oxygen saturation may then be expressed as: 



SO, 



k^ii (9c) 



Other variations of the foregoing methods leading to the equations (9b) and (9c) may also be 
used as long as they are designed to eliminate system parameters and to ultimately express 
[Hb] [HbO], and/or oxygen saturation in terms of known or measurable system variables or 
10 parameters such as, e.g., optical densities, extinction coefficients, and/or other geometry- 

dependent parameters. 

It is further appreciated that the foregoing method of the present invention allows the 
wave sources to irradiate multiple sets of electromagnetic waves having different wave 
characteristics through various different embodiments. The simplest arrangement may be to 
15 provide two wave sources (such as SI and S2), where each source is designated to irradiate 
electromagnetic waves having different wavelengths, phase angles, amplitudes, and/or 
harmonics. Alternatively, each wave source may also be arranged to irradiate substantially 
identical signal waves which are, however, superimposed on different carrier waves. In yet 
another alternative, a single or each wave source may be arranged to irradiate multiple sets of 
20 electromagnetic waves intermittently, sequentially or simultaneously as long as different sets 
of electromagnetic waves can be identifiable by one or more wave detectors. Similar 
arrangements may also be applied to the wave detectors as well. For example, two wave 
detectors (Dl and D2) may be provided where each detector is designated to detect only a 
single set of electromagnetic waves. Alternatively, a single or each wave detector may be 
25 arranged to detect multiple sets of electromagnetic waves with different wave characteristics 
on an intermittent, sequential or simultaneous mode. Because the foregoing methods of the 
present invention allow these various arrangements, the foregoing methods can be readily 
incorporated to any conventional spectroscopic techniques including, e.g., continuous wave 
spectroscopy, time resolved spectroscopy, and phase modulation spectroscopy. 
30 " In another aspect of the invention, an over-determined numerical method is provided 
to solve the modified Beer-Lambert equation and/or the photon diffusion equation applied to 
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an optical system including a source module and a detector module, where at least one of the 
source module and detector module may be arranged to irradiate or detect more than two sets 
of electromagnetic waves. By arranging an optical system to provide more equations than 
the number of system variables of interest, resulting extra equations may be utilized for other 
5 purposes, e.g., (i) to enhance the accuracy of estimated values of system variables (e.g., 
chromophore concentrations or their ratios), (ii) to determine system parameters (e.g., «a m ,» 
«(3 n ," w y," w Bmn," "W "5," "Ei," "a" or other parameters such as absorption and 
scattering coefficients of the medium and/or chromophores) or (iii) to provide correlations of 
the medium- and/or geometry-dependent parameters of the equations (1) or (3b) with the 
1 0 system variable(s) and/or other system parameters. 

In the first embodiment, the extra equations may be used to obtain multiple sets of the 
chromophore concentrations (and/or their ratios). It is expected that discrepancies may exist, 
at least to some extent, among the estimated values of the concentrations (and/or their ratios). 
Such discrepancies may be attributed to inherent idiosyncracies of each pair of the wave 
15 source and the wave detector. Alternatively, the discrepancies may also arise from a non- 
homogeneous medium, i.e., the medium having inherent variations in optical properties in 
different portions thereof. One way of taking advantage of different concentrations of the 
chromophores (and/or their ratios) may be to average such values to obtain an arithmetic, 
geometric or logarithmic average having reduced random or systematic errors and greater 
20 accuracy. Alternatively, each measured value may be weighted by an appropriate weight 
function which may account for, e.g., geometric configuration of the wave source-detector 
assembly. 

In the second embodiment, correlations of the medium- and/or geometry-dependent 
parameters of the equations (1) or (3b) with, e.g., the chromophores concentrations (or their 
25 ratios) may be obtained from those extra equations. For example, when G (i.e., the ratio of 
F \ to F ^ ) is approximated as a polynomial of oxygen saturation according to the equation 
(10), each of the coefficients of the polynomial may be assigned an initial value, and then 
improved by iterative techniques employing a conventional numerical fitting method. In 
addition, the extra equations may also be used to find the correction functions between 
30 approximated values of oxygen saturation, [Hb], and/or [HbO] and the true values thereof. 
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Furthermore, the extra equations may also be used to estimate system parameters 
(e.g.,«a m ,» «P„,» "Y,» "B^," "5," «a," and/or other system parameters 

such as absorption coefficients and/or scattering coefficients of the medium and/or 
chromophores). For example, a forward numerical scheme may be used to estimate 
absorption and reduced scattering coefficients of the physiological medium and/or 
chromophores included therein. As described hereinabove, migration of photons and 
propagation of electromagnetic waves in the medium can be described by the diffusion or 
transport equation. Assuming that the medium is semi-infinite and homogeneous, following 
equation may describe an intensity of electromagnetic waves detected by a j-th detector: 
I \ (13) 

where Si generally denotes a source coupling parameter accounting for, e.g., characteristics 
of an i-th wave source such as power and configuration thereof, mode of optical coupling 
between the i-th wave source and medium, and/or coupling loss therebetween, and Dj is a 
detector coupling factor generally accounting for characteristics of a j-th wave detector, 
mode of optical coupling between the j-th wave detector and medium, and the associated 
coupling loss therebetween. 

A symbol represents a forward numerical model which simulates measurement 

for a given pair of the wave source and wave detector. Parameters «p a " and «m>" 
represent, respectively, an absorption coefficient and reduced scattering coefficient. When 
the optical system includes, e.g., a total number of N s wave sources and N D wave detectors, 
equation (13) can be expressed in a matrix form as follows: 

(14) 



25 











A, • 









Each side of the equation (14) is divided by the first column of each matrix: 



30 
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Each row of the matrices 
matrix: 



ices A and B of the equation (15) is then divided by the first row of each 



A = 



I 



N S .N D 



I 



N s ,l 



^(r Ns ,r„/i a ,^) K r i' r i'^^s) . 



= B 



(16) 



As manifest in equation (16), both of the matrices A and B are functions of the absorption 
and reduced scattering coefficients and do not depend on the source- and detector-coupling 
parameters such as S, and D, Accordingly, by minimizing the difference between two 
matrices A and B (i.e., || A - B || ), the best estimates of the absorption coefficient and the 
reduced scattering coefficient may be numerically obtained by conventional curve-fitting 
methods. After estimating the absorption and reduced scattering coefficients, [Hb], [HbO], 
and oxygen saturation may be obtained by the following set of formulae: 



[Hb} = 



[HbO] = 



*HbO< 
£ Hb e HbO 



C J 1 n** _ £ ** ll k 



(17a) 
(17b) 



[HbO] 



s ° 2 = [m]+[Hbo] = fat - (-itf - £ >>) 



(17c) 



It is noted that the foregoing over-determined method may be applied to the optical 
systems with at least two wave sources and three wave detectors, at least three wave sources 
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and W „ wave detectors, or three wave sources and three wave detectors. AHernattvely, the 
over-determined method may equally he applied to the optical systems where a single or each 
wave source or detector has the capability of irradiating or detecting multiple sets of 
electromagnetic waves, respectively. 

It u appreciated that the foregoing over-determined method may be incorporated by 
any conventional numerical scheme, For example, a forward model, a backward mode, or a 
hybrid model may be applied to determine, e.g., an extinction coefficient, an absorption 
coefficient or a scattering coefficient of the physiological medium (or the chromophores 
included therein). Such models may also be applied to estimate the concentrations of the 
chromophores (and/or ratios .hereof). It is uoted, however, that the results obtained by such 
numerical models genmlly include inherent errors associated therewith. Such errors may be 
minimized by employing numerical models with the error terms of the second or htgher 
order However, such models may have a major drawback of requiring rigorous numencal 
computation, Accordingly, the accuracy and efficiency of each numerical model must be 
considered and compromised in selecting an appropriate one. 

In yet another aspect of the invention, an optical system is provided to solve a set of 
wave equations and to determine absolute values of the concentrations of the chromophores 
(and/or ratios thereof) contained or suspended in a physiologica. medium. An exemplary 
optical system may include a body, a source module including a, .east one wave source, a 
detector module having a, .east one wave detector, and a processing module. The source 
module is supported by the body, and is generally arranged to optically couple wtth me 
physiological medium and to irradiate into the medium at leas, two sets of e.ectromagnetrc 
waves having different wave characteristics. The detector module is also supported by the 
body and is arranged to optically couple with the medium and to detect electromagnettc 
waves transmitted through the medium. The processing module is arrange* to operattvely 
couple with the detector module, to solve a set of multiple wave equation, and to obtam the 
absolute values of the concentrations of chromophores and/or ratios thereof. 

In general, the processing module includes an algorithm which is arranged to solve 
the foregoing equations (1) or (3b) or their modified versions. For example, one or more of 
, the foregoing methods may be provided thereto as hardware or software or be implemented 
,o a microprocessor so that the absolute vataes of the chromophore concentrations (and/or 
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their ratios) may be calculated from, e.g., the intensity of electromagnetic waves irradiated by 
the source module (or its wave source(s)), intensity of electromagnetic waves detected by the 
detector module (or its wave detector(s)), and at least one system parameter which may 
account for an optical interaction and/or interference between the electromagnetic waves and 
5 the medium. The algorithm of the processing module may include one or more functions 
which correlates the medium- and/or geometry-dependent term(s) of the foregoing wave 
equations as a function of, e.g., the chromophore concentrations (or their ratios). The 
algorithm of the processing module may also be capable of executing the over-determined 
methods described hereinabove. In addition, the processing module and algorithm thereof 
10 may be arranged to operate on the TRS, PMS or CWS mode. 

The source module may include at least one wave source and the detector module 
may include at least two wave detectors. Alternatively, the source module may include at 
least two wave sources while the detector module may include at least one wave detector. It 
is preferred, however, that the source and detector modules include, respectively, at least two 
1 5 wave sources and at least two wave detectors. 

As described hereinabove, the foregoing methods of the present invention are rather 
insensitive to actual configuration of wave sources and detectors. Accordingly, the optical 
system of the present invention may include any number of wave sources and/or detectors 
arranged in any arbitrary configuration. However, a few source-detector configurations may be 
20 preferred so as to obtain the absolute values of the chromophore concentrations (and/or their 
ratios) with better accuracy, reliability, and reproducibility. 

In the first preferred embodiment, multiple wave sources and/or detectors may be 
arranged so that near-distances between each wave source and detector pair are substantially 
identical. For example, for the source module including two (i.e., a first and second) wave 
25 sources and the detector module including two (i.e., a first and second) wave detectors, a 

distance (e.g., a near-distance) between the first wave source and the first wave detector may 
be arranged to be substantially similar to that (e.g., another near-distance) between the 
second wave source and the second wave detector. In addition, another distance (e.g., a far- 
distance) between the first wave source and the second wave detector may also be arranged 
30 to be substantially similar to that (e.g., another far-distance) between the second wave source 
and the first wave detector. It is appreciated that this embodiment is not necessarily required 
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for every single pair of the wave source and detector. For example, when the source module 
has M wave sources and the detector module has N wave detectors (M and N are tntegers 
-eater than 1), a, least mo of M wave sources and two of N wave detectors may be arranged 
so mat a distartce between an M,-.h wave source and an N,-th detector is substantially stmilar 
to that between an M> wave source and an N r th wave detector, and that a dtstance 
between the M,-m wave source and the N,-th wave detector is substantially similar to tha, 
between the M 2 -th wave source and the N,-th wave detector, whore M, and M, are both 
integers between 1 and M, and N, and N : are both integers between 1 and N. 

In the second preferred embodiment, wave detectors may be arranged substantially 
linearly along a straight line. As will be demonstrated in Examples, this arrangement has 
produced the concentrations of oxygenated and deoxygenated hemoglobins and oxygen 
saturation with great accuracies. It is noted tha. not all wave detectors have to be arranged 
linearly. For example, a substantial portion of the wave detectors may only have to be 
disposed substantially along the straight line. In this embodiment, it is also preferred that at 
leas, one wave source be disposed on one side of the straight line while at least one other 
wave source be disposed on the other side of the straight line. 

In operation, a source module with at least one wave source and a detector module 
having at leas, one wave detector are provided to an active surface of an optical probe whtch 
is operatives connected with a main body of an optical system. Alternatively, the wave 
source and/or detector modules may be disposed a. the main body and optical fibers may be 
provided to connect the source and detector modules to openings provided on the active 
smface of the optical probe. Any conventional wave sources and detectors may be used for 
such optical prove. It is preferred, however, tha. the wave source is capable of trradtating 
electromagnetic waves in .he near-infrared range, e.g., between 500 and 1,200 nm or, in 
particular, between 600 and 900 nm, and that the wave detector has appropriate sensitivity to 
such electromagnetic waves. The optical probe is placed on a physiological medium, w«h » 
active surface disposed on me medium to form an optical coupling therebetween. The source 
module is activated so mat a. leas, two sets of electromagnetic waves having different wave 
characteristics are irradiated into the medium. The detector module men picks up different 
, sets of electromagnetic waves irradiated by the wave source, propagated through the 

medium, and directed toward the wave de.ec.or. The wave detector generates eleotnc signals 
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which are received hy the processing module of the main body of .he optical system. Based 
on the optica, densities and a. leas, one system parameter snch as extinction coefficients of 
the chromophores, me processing module computes the absolute values of the oxygenated 
hemoglobin concentration, the deoxygenated hemoglobin concentration or the oxygen 
5 saturation. 

B is noted mat .he optical system according to me present invention may tnclude an 
equation solving module which is operationally separate from the processing module. Such 
an equation solving module may include variety of numerical models designed perform one 
or more of the foregoing methods of the present invention. 

0 Although the foregoing disclosure has been directed toward oblaining the absolute 

values of the concentrations of oxygenated and deoxygenated hemoglobins (and/or then 
mtios), the foregoing optical systems and methods therefor may be applicable to obtam 
absolute values of other substances in the medium or properties thereof. For example, 
similar or identical systems and methods may be used to determine .he absolute values of 

15 concentrations (or .heir ratios) of other chromophores such as lipids, cytochromes, water and 
the like. Depending upon the absorption or scattering characteristics, the wavelengths of the 
electromagnetic waves may be adjusted for better resolution. In addition, chemtcal 
compositions may be added to the medium to enhance optica, interaction or interference of 
chromophores in the medium or to convert an non-chromatic substance of the medium tnto a 

20 chromophore. 

As described hereinabove, the foregoing optical systems and methods of the present 
invention are preferred to be incorporated to the continuous wave spectroscopic technology. 
However, snch systems and methods may readily be incorporated into the time-reso.ved and 
phase-modulation spectroscopic technologies as well. 

25 The optical systems and methods according to the present invention find a variety of 

medical applications. As described hereinabove, such optical systems and methods may be 
applied to measure the absolute values of concentrations of oxygenated and deoxygenated 
hemoglobin and/or their ratio. Such optical systems will be beneficial in non-invastvely 
diagnosing ischemic conditions attd/or locating ischemia in various organs and tissues such 

30 as e g a brain (stroke), heart (ischemia) or other physiological abnormalities ongtnattng 
from or characterized by abnormally low concentration of oxyhemoglobin. In addition, 
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presence of cancerous tnmors in varions internal organs, breasts, and skins may be eaady 
detected as well. Such optica, systems and methods may further be applied to cells disposed 
in epidermis, corium, aod organs such as a tang, .iver, and kidney. Such optica, systems and 
methods may also be applied to diagnose vascular occlusion durhtg or after surgrcal 
procedures including transplantation of tissues, skins, and organs, e.g., heart, tang, hver, and 
kidney. 

Other medical applications include tissue analysis of regions and/or of an location 
over time FIG. 3 shows an embodiment of an optical system of the invention that includes 
multiple probes. A computer system 1 is in communication with multiple probes 201. 
Although three probes are shown, embodiments can have more or fewer probes. 

Computer system 1 directs the electromagnetic waves that are generated by each 
probe and receives data on the electromagnetic waves that are detected in for use in 
calculation, In one embodiment, the probes have two wave sources and four wave detectors 
as shown in FIG. 4. 

The multiple probes can be perform simultaneous measurements at different locatrons 
on a physiological medium for self calibration. Additionally, the measurements from the 
multiple probes can be utilized to generate regional information about (he physiologrcal 
medium The regional information can include 2D mapping of me surface of the medrum. 

One or more probes can be utilized to monitor changes in the physiological medtum 
over time. This can be useful for medical applications such as monitoring tissue viabthty. 
For example, in one embodiment, me probes are detachable and flexible, altawing them to be 
p.aced on .he skin of a patient after surgery, such as plastic surgery. The oxygen saturation 
in tissue can be utilized as an indication of tissue viability. 

The multiple probes can be utilized to calculate the tissue viability at different 
; locations. Also, the probes can be left on the skin to allow monitoring changes in tissue 
viability over time. 

Following examples describe simulation and experimental results obtained by the 
optical systems and methods thereof according to the present invention. All simulation and 
experimenta. results indicate that the optical systems and methods thereof provide accurate 
0 predictions of the concentrations of the hemoglobins as well as the oxygen saturation. 
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F.X AMPLE 1 

The diffusion equation (3b) was numerically solved for optica, proves with multiple 
wave sources and detectors arranged in various configuration, The equations were apphed 
,„ a sample physio.ogica. medium such as a semi-infinite, homogeneous diffuse medtum 
with different backgmund opfieal properties. Diffuse reflectances were ca.cula.ed accordtng 
to an imaging source approach disc.osed in an arfic.e emitted, "Boundary conditions for .he 
diffusion equation in radiative transfer" by R.C. HaskeU e, a,, and pub.ished in Journal of 
Optica. SocietyofAmericaOI, p. 2727-2741, 1994). Values of G (the ratio of to f ) 
were estimated a, different levefs of oxygen saturation and fitted as a polynomtal .hereof. 

FIGa 5A to 5C are plots of G-s (ratios of F-a simulated a. different wavelengtits) as 
a function of oxygen sahuation As shown in me figures, all simulation results showed 
diatinc. one-to-one correlations between G and oxygen saturation, and ma. the dependence of 
G on oxygen saturation was substantial., insensitive to ttte source-detector configuratton as 
wel. as me background optica, properties. Conventional curve-fitting methods, such as the 
.east-squares memod, were applied to numerical., estimafe the coefficients of equation (.0) 
<■*.**,*....> For example, in a aysfem with background reduced acattenng 
coefficient of 10 em" and me torn, hemoglobin concentiation of 10 J mol/mer and s,mu,a,ed 
with elecfromagnetic waves having wavelengtits of 780 nm and 830 nm, fo.lowing equation 
approximated the relation between G and oxygen saturation: 

(18) 

n = — - 0 728 + 0.399 SO, * 0.064 • SO 1 , + 0.067 • SO, 
F h 

FYAMPLE 2 

Further simulations were performed in a system with the background scattering of 7 
cm" and the total hemoglobin concentration of 2x10^ mol/liter. In the simulations, S0 2 was 
varied from 0 to 100%. FIG. 6 is a plot of calculated oxygen saturation contrasted agaxnst 
true oxygen saturation. Although the background properties used to find the correlate 
between G and oxygen saturation were quite different, the estimated oxygen concentration 
was accurate with a systematic error of about a few percent. 



) 



EXAMPLE 3 
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An exemplary optica, system was prepared and hemog.obin concentrations and the 
oxygensaturation were monitored before and after occlusion of arteries of a. extremfty m a 
human subject. FIG. 7 is a schematic diagram of a sample optica, system including .wo 
wave sources and two wave detectors arranged in a Hnear fashion acceding to the present 
invention. Two wave detectors were linearly disposed a,on g a straight Une and separaed by 
6 nun. Two wave sources were disposed outside of each wave detector so that the left wave 
source was disposed a, the .eft side of the .eft wave detector at a distance of 9 nun, and the 
ngh, wave source disposed a. the right side of the right wave detector a, a drstance of 9 mm. 
According each pair of me wave sources and detectors had identica, near-distances sard 
far-distances. The wave sources had an outer diameter of 2 mm and each was designed to 
toad-ate the electromagnetic waves havmg two different wavelengths. For example, aser 
diodes (model numbers HL6738MG HL8325G, Thorlabs, Ine), were used to rrradrate 
efcctromagnetic waves having wavc.engths 690 nm and 830 nm, respective*. A photo- 
detecto, (model number OPT202 by Ban-Brown) was used as the wave detector. 

A euffwas placed at the upper arm and the optica, probe was disposed a. the fore 
m After the subject stabilized, cuff pressure was increased to about .60 mmHg in about 
35 seconds, he.d a, the same level for about 40 seconds, and men re.eased to atmospheric 
.eve.. Concentrations of the tota, hemoglobin, oxygenated hemoglobin, and deoxygenated 
hemoglobin were monitored, along with the oxygen saturation. 

FIG 8 is a time-course plot of total hemoglobin (HbT) concentration, oxygenated 
hemog.obin (HbO) concenbation, and deoxygenated hemoglobin (Hb) concentratton 
according to the present invention, and FIG. 9 is a time-oourse p.o, of oxygen sa.ura.ton 
according to the present invention. As shown in the ngrrres, hemogmbin concensus and 
oxygen saturation decreased sharply during the initial phase of occlusion, fo.lowed by a 
gradual decrease thereof. After the release, concentrations and oxygen saturation showed 
Id increase. These resuhs demonstrated mat the optical systems and methods accordmg 
th e present invention provided accurate predictions of the hemoglobin concentrations as weft 
as the oxygen saturation. These results also showed mat me optical systems possessed 
proper temporal response characteristics. 
, It is to be understood mat, while various embodiments of une invenfon has been 

described in conjunction with the detailed description thereof, the foregoing is intended only 
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,0 fflustrate and no, «o limit the scope of d. invention, whieh is defined by the scope of the 
appended claims. Other embodiments, aspects, advantages, and modifications are wilhan me 
scope of the following claims. 



PHOTP003 



32 



